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I present a simple discussion of the masses and mixings of the light pseudoscalar and vector 
mesons based on a l 'qq" description of the effective field theory. The analysis includes 7/(958) from 
the beginning, and is largely concerned with structural questions. While the final results are mostly 
known, the method gives insight into the general form of the meson mass matrices and the different 
character of the mass splittings and mixings in the pseudoscalar and vector multiplets, and provides 
a coherent overall picture 
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The theoretical analysis of the masses and mixings of the light pseudoscalar and vector mesons has a long history, 
from the original Gell-Mann-Okubo mass formula through modern analyses based on chiral perturbation theory 
as summarized, for example, in The masses of the light pseudoscalar mesons are now discussed in terms of the 

■ breaking of the chiral SU(3) L x 811(3)^ symmetry that exists in QCD for massless quarks, with the octet mesons 
appearing initially as the eight Goldstone bosons for the broken symmetry. The SU(3) singlet state |?7i) ~ \rj'} and 
the vector mesons are treated separately. 

A perturbative treatment of the effect of quark masses on the pseudoscalar mesons, introduced through the quark 
mass matrix 
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I. INTRODUCTION 



"pH , leads to the prediction that HQ 
Oh! 

m; g = -(2m s + m)B , m 2 m = m 2 + -(m s + 2m)B . (1-2) 



m, = IB^rh, m 2 K — (m s + m)Bo, 



Here fh = {ra u + m,d)/2 is the average light-quark mass, m 2 K = {m 2 K - + m /fo)/2, and ra^ and m ril are the masses 
of the unmixed SU(3) isospin-0 octet and singlet states \r)%) and \r)i), Bq is an unknown dynamical matrix element, 
and mo is the original mass of the 771, also unknown. Finally, the original states \r]s) and \rji) couple through an off 
C$ I diagonal matrix element 

2 2y2 . 

m 8,i = —(m s -m)B Q . (1.3) 

to produce the observed mass eigenstates |ry) = 77(547) and |r/) = 77(958) 0]. 



The perturbative predictions are reasonably successful. Equation (1.2) contains the Gell-Mann-Okubo mass relation 



m 2 Vs = I (4M| - ml) , (1.4) 
which is valid to about 20% when the 77s is identified with the 77. The 77-77' mixing angle defined by 

I?/') = sinflp 1 77s) + cosflp \r)i) , \rj) = cos6 P \r] S ) - sind P \t]i) , (1.5) 
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is predicted to be —10° S, at the edge of the range —13° to 22° obtained in various fits to decay data, for example, 

The situation is different for the vector mesons, where the chiral symmetry arguments do not apply. However, 
symmetry ideas from SU(6) and quark-model and similar methods again lead to successful mass relations among the 
vector mesons, while dynamical ideas related to annihilation processes have been invoked to explain the different 
mixing patterns fll],f§ . 

In the present paper, I present a simple structural analysis of the mass and mixing problems for the nonets of 
pseudoscalar and vector mesons. Most of the ideas appear in the literature as applied to specific problems, and most 
of the results are also familiar. However, by introducing some quar k-mass effects which modify previous approaches 
and lead naturally to mass matrices of the most general type |12| , I find a consistent overall picture in which, for 
example, the different mass and mixing patterns for the pseudoscalar and vector mesons follow from physical input, 
and the ratios of dynamical matrix elements in the two multiplets give information on the mixings. I first treat 
the pseudoscalar nonet, with the rf included from the outset. I then treat the vector nonet, and comment on the 
connections and differences, and conclude with a few overall remarks. 



II. THE PSEUDOSCALARS 



A. Structure of the mass matrix to 0(m 9 ) 

In this section, I discuss a generalized version of the mass matrix which allows actual fits to the masses of the 
pseudoscalar mesons. It will be useful to start in a "quark model" basis of effective meson fields <fiij{ x ) an d meson 
states \4>ij), where i, j label the flavor structure in a 3 x 3 description of a flavor SU(3). The general basis includes 
both the SU(3) octet and singlet states. The mass operator for the pseudoscalar mesons is then of the form 

M o P = E <^f M h'^r 

I will work throughout with the matrix M 2 , and will first consider the neutral sector. 

The combinations uu, dd, and ss for the flavor labels i, j are completely equivalent for massless quarks in QCD, 
so, with a suitable choice of phases, the initial mass matrix for the neutral states must have the same entry in every 
location corresponding to an SU(3) singlet configuration, 

^ / rriQ rriQ m\ \ / mi mi mi 

M op = o m o m o m o \ = B o ™i mi ™i I ■ (2.2) 



3 



,2 ™2 



mi mi mi 

The factor Bq in the second form has been extracted for convenience in writing later equations. Mq P has eigenvalues 
zero in the SU(3) octet states 1 7r) and |t7s), or any combination thereof, and the eigenvalue m§ = SB^mi for the 
SU(3) singlet |?7i). The pion and 77s are therefore massless, as expected for Goldstone bosons. The term m\ contains 
contributions from the a xial anomaly [Q and from annihilation diagrams such as that in Fig. [l](a). Both lead to the 



singlet structure in Eq. (2.2) when written in the qq basis. 

The structure changes in the presence of quark masses, both through the addition of diagonal contributions pro- 
portional to the quark masses, and through the effect of those masses on the annihilation terms. The effect of the 
latter is neglected in all discussions of which I am aware. The modified mass matrix has an obvious structure 

(mi + 2m u + 2am u mi + a(m u + 1714) mi + a(m u + m s ) \ 
mi + a(m c i + m u ) mi + 2md + 2amd mi+a(md + m s ) (2-3) 
mi + a(m s + m u ) mi + a(m s + rrid) mi + 2m s + 2am s I 

in the qq representation, where the terms proportional to a arise from mass insertions in the annihilation diagrams.^ 
I expect a to be small, with |a|m s <^ mi. A dispersion relation argument shows that the annihilation contribution 



x The basic structure of Mp was apparently first suggested by Isgur juj, but without the a terms. Isgur identified what are 
here the Bomi terms with the annihilation diagrams, treated the diagonal terms phenomenologically as effective quark masses, 
and used the structure to explain the difference between the rj, 7/ and ui, (f> mixings. The general form of the mass matrix was 
given by Morpurgo 112] , but without the physical identification of the various contributions. 
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to mi is positive in the symmetrical limit. Its magnitude is decreased by mass insertions in the loop, so a should be 
negative and small.^] The light quark masses are also small, m u , md <C m s . 

A transformation to the SU(3) basis gives a considerably less transparent structure. The resulting matrix can be 
simplified by dropping a diagonal term proportional to am and off-diagonal terms which involve m u — md- Both are 
small, with \am\, \m u — rrid\ <C m s , mi. The off-diagonal terms induce a small mixing of the ir° with the rj and rf , 
but first contribute to the physical masses in second order . The reduced matrix is 



2m 

Mg=B Q \ 2m + 4m s /3 -V2(2 + 3a)m s /3 

-V2(2 + 3a)m s /3 2m + 3mi + 2(1 + 3a)m s /3 



(2.4) 



where m = (m u + md)/2 and m s = m s — m. It is convenient to rewrite this as 



Ml = mil 





T M 2 



(2.5) 



where 1 is the 3x3 unit matrix, is the 1x2 null matrix, and Mp is the 2x2 matrix 



Ml 



4m s /3 



-V2(2 + 3a)m s /3 



-y/2(2 + 3a)m s /3 3mi + 2(1 + 3a)m s /3 



(2.6) 



The matrix m^l contributes a term m\ to each of the eigenvalues of Mp. I will drop this term and work with the 
2x2 matrix Mp which has eigenvalues m 2 = m^ — m\ and m^, 
Mp has the form of the general ijs, r\\ mass matrix 0,0, 



M z , = 



Mis M x 2 8 
M? 8 Ml 



(2.7) 



with Mf 8 real after a phase has been absorbed. In particular, both involve three independent parameters. I find it 
striking that the explicit inclusion of the annihilation terms already leads to this structure at first order in m s . 

The kaon and charged pion masses have the usual simple structure which does not involve annihilation contributions, 
with 



m 2 K+ 



B (m s + m u ), m 2 K o = B Q (m s + m d ), m\ ± = B (m u + m d ) = 2B m. 



The average kaon mass is m 2 K = (m 2 K+ + m^._)/2 = Bo(m u + m), so 



ml = m 2 K = B m s 



(2.8) 
(2.9) 



B. Physical masses and fits at 0(m q ) 



The masses of the physical particles rj and Vj are determined by the eigenvalues of Afp, 

m 2 = P m i + 2(1 + a)m s ] — t^B® [(^ mi + ^(1 + a)m s ) 2 — 8m s (2mi — a 2 m s )] 



1/2 



ih 2 , 



4 

^B m s 



, 2/ 3 \ m s 
1 - - 1 + -a — 
9 V 2 I mi 



^B [3m! + 2(1 + a)m s ] + ^B [(3m 1 + 2(1 + a)m s ) 2 - 8m s (2m 1 - a 2 m s )] 



1/2 



3-B mi + -B m s 



4 / 3 \ m s 
l + 3o+-l + -o — 
9 V 2 J mi 



(2-10) 



(2.11) 



2 The rj' acquires mass through the mi terms without any quark mass insertions. As a result, the off-diagonal contributions 
proportional to m a can be complex, with Mp still Hermitian. The imaginary part of a should be small, and can be absorbed 
to first order in a change in the phase of |rji). I will therefore take a as real. 
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where ml, = ml, — to 2 and to 2 = to 2 — m\. The leading terms give the usual first-order perturbative result except 
or the independent appearance of a in the expressions for mt, and the rj, rj' mixing angle. 
If a is small, we obtain the inverse relations 

1 / 1 



B m s 



3^2 



B mi pa - | in 



2< 



(2.12) 



for Bom s and B^rnx. The first becomes the Gell-Mann-Okubo relation if Bom s is identified with rh 2 K as in Eq. ( p.9| ). 
These relations give the numerical estimates Bom s ps 0.211 GeV 2 , £>o m i ~ 0.252 GeV 2 , and m s /mi ~ 0.8 73 fo r the 
physical masses to 2 , = 0.898 GeV 2 , to 2 = 0.281 GeV 2 . The expansion parameter 2m s /9mi pa 0.19 in Eqs. ( 2.10| ) and 
( [2.11 ) is therefore small, but not so small that the correction terms can be ignored. Since the "first order" fit to the 



observed masses automatically incorporates these corrections and other diagonal m g terms in the fitted values of m\ 
and M s , it seems preferable to use the general form of the eigenvalues or Mp in physical analyses. 
The eigenvalues of the general mass matrix Mp are determined by the conditions 



m 2 , + m 2 = TrM| = B [3mi + 2(1 + a)m s ] 



m\,m\ 



det Mp = B^m x m s - 2a z m 2 ). 



(2.13) 
(2.14) 



The a 2 m 2 term in the last expression is presumably small, and will be dropped in the immediately following equations 
for simplicity.^] 

Eliminating m s between the two equations and solving for toi gives the relation 



Bomx = i (m 2 , + to 2 ) + i (to 2 , + m 2 ) 2 - 6to 2 ,to 2 (1 + a) 



1/2 



(2.15) 



The argument of the radical is negative for the physical masses, and there is no real solution for toi, unless 



a < 



to ?j , + m„J — 6m^/TO^ l&rh^m^ = —0.084 



(2.16) 



This gives a useful constraint on o ur ge neral parametrization of Mp and requires a to be ne gativ e, as expected. 
Once a value of a is given, Eq. (2.15) determines Borrii, and Bom s then follows from Eq. (2.14), 



B m s = to 2 ,to 2 /4S toi 



(2.17) 



with possible corrections of relative order o 2 to s /2toi. 

We n eed one further condition to determine a. This is given within the rji, rj% sector by the mixing angle 9p defined 
inEq. (0), 



tanflp = -2V2(2 + 3a) to s /{9toi - 2 (1 - 3a) m s 



(3mi + 2m s + 2am s ) — 8 (2toi — a 2 m s ) to< 



1/2 



(2.18) 



This angle is not well determined experimentally. Different analyses give values of —13° to —22° depending on the 
data used p|,|6|-|Io|| , and the theoretical connections assumed to hold between decay matrix elements and the meson 
masses and mixings. 

A different approach is to take the value of Bom s w m 2 K as known. The parameters toi and a are then determined 
by the rj and rj 1 masses, and Bp can be predicted. I will adopt this approach, but will allow for some difference between 
the effective values of Bom s in the kaon and 77, 77' systems by introducing a multiplicative correction factor x, with 
BQm s = xrh 2 ^. This difference can arise from higher order effects in the chiral expansion, loop corrections, or the 
change in matrix elements expected in the change from a system with one heavy quark to systems which can have 
two. All of these corrections can be encompassed in our general parametrization of Mp; their separation depends on 



It w ill actually be small for the final range of parameters discussed below. Including this term pushes the limit on a in Eq. 
(2.16) farther from zero. 
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dynamical calculation beyond the scope of this paper. The net effect gives x ^ 1, with the expectation that x — 1 is 
on the order of the fractional violation of the Gell-Mann-Okubo relation. 

The results obtained in our approach are summarized in Figs. || and [3] which show the fitted value of a and the 
predicted mixing angle 9p as functions of x. It appears from the results on mixing angles that x is in the range 
1.0 < x < 1.2, a reasonable result, with corresponding values of a in the range — 0.2 4 < a < —0.1. The parameter o 
is therefore small, as expected. An argument given later in connection with Eq. (3/7) suggests that, when the vector 
mesons are included in the analysis, values x w 1.05, a rs —0.18, and 6p w —15° are favored. This predicted mixing 
is in the center of the ranges found in various analyses [BMlO) ■ 



III. THE VECTORS 



The situation for the vector mesons p, to, and 4> is quite different. In the qq basis, the initial mass matrix has 



the form of the pseudoscalar mass matrix in Eq. (2.3) with the addition of a diagonal term 5jmyl which gives a 
common mass to all the vector states. No such term was not allowed in the pseudoscalar sector by the requirement 
that the pion and r)% appear initially as Goldston e bo sons of t he broken SUl(3) x SUr(3) symmetry. Furthermore, the 
elements analogous to the terms B$m\ in Eqs. ( |2.2| ) or (2.3) have no anomalous component, so obtain contributions 
only from annihilation diagrams. The leading vector annihilation diagram involves three-gluon exchange as shown 
in Fig. ^(b). The corresponding matrix element, which I will denote by B' m 2 , is expected to be much smaller than 
the two-gluon contribution to the pseudoscalar masses from Fig. 0(a), but still positive with 7772 171%. The relation 
mi ~ m s from the pseudoscalar sector then implies that 7772 <C m s . The effect of mass insertions on the already-small 
m2 terms should also be small, with |a'| <C 1. Dropping those corrections, we obtain the simplified vector mass matrix 



My 



my + 777*2 + 2m u m 2 m 2 

B' n I 777 2 my + m 2 + 2m d m 2 

777 2 m 2 my + 777 2 + 2m s 



(3.1) 



The matrix elements Bom s and B' m s differ in the dynamical quark model only by the effects of spin-dependent 
interactions on the wave functions and are expected to be similar in magnitude, Bq ~ B' . 

In contrast to Mp, which was dominated by the flavor-singlet structure of the mi terms in Eq. (2.2) with 2m s < 3mi 
appearing as a perturbation, My is dominated by the the strange-quark mass term. There is consequently no reason 
to transform to the SU(3) basis used to pick out the singlet structure. The qq basis is natural. It is useful, however, 
to transform to a flavor SU(2) basis \cj>i) — (\4> U u) — \4>dd)) /v2, |</>o) = (|0«a) + \4>ddl) /v2 for the light-quark sector. 
This isolates the isovector component \<f>\) = \p), and brings the mass matrix to the form 



My = B' (m v + 2m) 1 





T My 



where I have dropped off-diagonal terms in ?77 u — 777^ which contribute linearly to the p, to mixing 
quadratically to the masses. The matrix My is 



(3.2) 
but only 



My = B'q ( ^ ^™ 2 

° 1 V2TO2 2777,5 + 777 2 



(3.3) 



where m s = m s — m as before 
The mass of the p meson is 
the small quantity m 2 /m s , are 



The mass of the p meson is clearly given by m 2 p = B' (m v + 2tt7). The eigenvalues of My, taken to first order in 



777,., = 777,., - 771, 



= 2B' m 2 and 



= B' (2m s 



m 2 ) 



(3.4) 



These relations give B' (m v + 2m) = 0.592 GeV 2 , 2B^m 2 = 0.0206 GeV 2 , and 273^ = 0.437 GeV 2 . Thus, 
fn 2 jm s w 0.047, and 7772 <ti m s as expected. The approximate degeneracy of the p and u> masses follows from the 
smallness of 7772. This is expected in QCD as argued above. There is no constraint on the extra parameter a' in the 
general case since m 2 <C (m 2 + m|) 2 . 

The average K* mass squared m 2 K * is defined by 



L K , 



l K+ 



/2- 



B' m s 



(3.5) 



The combination of the two preceding equations gives the sum rule 
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s 2 x - 2 

m 4> ~ 



2m 2 



A"* 



(3.6) 



This is satisfied to 6% (0.437 GeV 2 « 0.412 GeV ). However, this sum rule is subject to some uncertainty theoretically 
because of different effect of the 0(m 2 ) contributions and loop corrections on the oj, 4> and K* systems. If the 

the scale factor x' is 1.06, in 



'A 



uncertainties are incorporated in a scaling of the K* contribution with B' m s 
the favored range for the scale factor x f or the pseudoscalar mesons. 

A comparison of Eqs. (2.12) and (3.4) leads to the rough estimate B /B' w 3m^/(2m? — to 2 ) ~ 0.96 while the 
ratio m 2 K /m 2 K , gives Bq/B' w 1.10. With either estimate, Bq/B' q 1 as expected. If I include the leading correction 
to TO 



from Eq. ( p. 10 ), I find the more precise relation 



Bn 



B' n x'm 2 



K 



3to 2 



A* 



2m 



-2 



2 3 \ m s 
1 - - 1 + -a) — 
9 I 2 J toi 



(3.7) 



This is satisfied for x — 1.05, a = —0.18. The value of x is in the range 1 < x < 1.2 consistent with analyses of to rj 
mixing, and gives the value Qp — —15°. The whole picture is therefore consistent. 

I will define the mixing angle 9y connecting the isospin zero state \4>o) with \(j> S s) so that 



\uj) = cosd v \<f> ) - sm9 v \<p sS ) 



sin9 v |0 O ) + cos6V \(f> s g) ■ 



The results above give 



0v 



TO 2 



V2r, 



TO 2 

2m s 



(3.8) 



(3.9) 



Numerically, 6y ~ m2/^/2m s w 1.95°, so there is very little mixing. The angle 6' v defined with respect to the SU(3) 
basis ^ is related to 9y by 



arctan(l 



35.3°, 



(3.10) 



where arctan(l/v / 2) = 35.3° is the angle for so-called "ideal" SU(3) mixing. The sum gives 8' v ~ 37.2°, in good 
agreement with experiment. However, I find little reason to use the SU(3) description. The inequality m s 3> to 2 picks 
out the u qq" basis states \4>ij) as giving the natural description of the problem, with the p split off in the isospin triplet 
substate. The physical masses then correspond yield a small mixing angle 9y ~ 1-9° without special assumptions. 

It is interesting to compare this result to that for the pseudoscalars. In the latter case, m s ~ Toi, so the SU(3) 
singlet is not clearly singled out, and it is less certain what basis is to be preferred. Mp is given in the \(f>o), \<f> s g) 
basis by 



2mi 



\/2(toi + arris) 



\/2(toi + am s ) m\ + 2(1 + a)m s 
The 0o, 0ss mixing angle 9' P defined by the relations 

|r/) = sin 9'p \4> ) + cos6»p \4> sS ) , |r/> = cos6»p |^ ) - sin6» P \(f> sg ) 



(3.11) 



(3.12) 



is given by 



tan#p = 2^2 (toi + am s ) / {—mi + 2 (1 + a) to 
+ (3toi + 2m s + 2am s ) 2 — 8 (2toi 
« \[2 (mi + am s ) / [mi + 2(1 + a)m s ] 



a 2 m s 



1/2 



(3.13) 
(3.14) 



This gives a rather large mixing angle in the range 47° to 33° for x from 0.9 to 1.2, a consequence of the similar sizes 
of toi and m s . The SU(3) basis is therefore somewhat better physically. In contrast, given the result m s ~ toi for 
the pseudoscalars, there is a clear physical argument that to 2 <ti mi ~ m s , so that the 4>ss basis is preferred for 
the vectors, with the mixing small corresponding to "ideal" mixing in the SU(3) basis. 
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IV. CONCLUSIONS 



My conclusions are simple: there is a substantial advantage in treating the nonets of pseudoscalar and vector meson 
together in analyzing the patterns of masses and mixings. It is important, in particular, to include the ?y'(958) in the 
analysis to get a consistent overall picture. Most of the results presented here are well-known, though not treated as a 
whole in other discussions. However, the complete analysis given here introduces a new feature through the through 
the inclusion of the leading quark mass correcti ons t o the annihilation diagrams in Fig. |l|. This leads to the general 
mass matrix for the ps eud oscala r me sons in Eq. (2.3), and a similar result for the vector mesons. These reduce to the 
simpler forms in Eqs. (2J5) and (3^2) if small off-diagonal terms proportional to m u — are neglected. The residual 
2x2 mass matrices have a completely general form. 

I show in the case of the pseudoscalar mesons that the physical r/ and rf masses can only be fitted if the annihilation 
parameter a in the mass matrix is nonzero and negative. The wide splitting of the tt, rj, and rj' masses, and fairly 
large mixing angle Op for the pseudoscalars, reflect the fact that the strange-quark mass insertion Bom s and the basic 
SU(3) singlet term B^mi are similar in magnitude. 

The situation is much simpler for the vector mesons. While the initial mass matrix has the same formal structure as 
that for the pseudoscalar mesons, the large anomalous contribution to the singlet component of the latter is missing, 
and the annihilation diagram is suppressed by the necessity of having at least three gluons in the intermediate state. 
As a result, the parameter mi analogous to m\ is expected to be small in QCD, with m,2 <C ni\ w m s , and the mass 
matrix simplifies. The smallness of leads to the near-degeneracy of the p and lo masses, and to the small mixing 
of the oj and 4> mesons in the natural qq basis. 

Finally, I show that the ratio of the matrix elements Bo and B' Q for the pseudoscalar and vector mesons is close to 
one, as expected. By calculating the ratio two ways, I find that overall consistency is obtained for a w —0.18 and scale 
factors x w x' « 1.05. The latter parametrize the higher order mass and loop corrections which differ for the K, K* 
and the rj, rf or oj, <f> sectors. The rj, rf mixing angle in the usual SU(3) basis is then predicted to be 6p = —15°. The 
oj, <fi mixing is small in the natural qq basis, Qy — 1.95°, so is nearly "ideal" in an SU(3) basis. 

I conclude that the whole picture is consistent. 



ACKNOWLEDGMENTS 



This work was supported in part by the U.S. Department of Energy under Grant No. DE-FG02-95ER40896. 



[1] M. Gell-Mann, phys. Rev. 125, 1067 (1962). 

[2] S. Okubo, Prog. Theor. Phys. 27, 949 (1962). 

[3] J. Gasser and H. Leutwyler, Phys. Rep. 87, 77 (1982). 

[4] J. F. Donoghue, E. Golowich, and B. R. Holstein, Dynamics of the Standard Model, (Cambridge University Press, Cam- 
bridge, England, 1992), Chaps. VII and XIII. 

[5] Particle Data Group, D. E. Grom et al, Eur. Phys. J. C 15, 1 (2000), Sec. 13.2. 

[6] J. F. Donoghue, B. R. Holstein, and Y.-C. R. Lin, Phys. Rev. Lett. 55, 2766 (1985). 

[7] F. J. Gilman and R. Kauffman, Phys. Rev. D 36, 2761 (1987); 37, 3348(E) (1988). 

[8] A. Bramon and M. D. Scadron, Phys. Lett. B 234, 346 (1990). 

[9] A. Bramon, R. Escribano, and M. D. Scadron, Eur. Phys. J. C 7, 271 (1999). 
[10] E. P. Venugopal and B. R. Holstein, Phys. Rev. D 57,4397 (1998). 
[11] N. Isgur, Phys. Rev. D 13, 122 (1976). 
[12] G. Morpurgo, Phys. Rev. D 41, 2865, (1990). 



7 




(a) 



(b) 



FIG. 1. Annihilation diagrams which contribute to the symmetric mass matrix Mq for (a), the pseudoscalar mesons, and 
(b), the vector mesons. 
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FIG. 2. Variation of the parameter a in the mass matrix Mp with the scale factor x in the relation Bom s = xm 2 K . Values 
of x in the range 1.05 < x < 1.2 are generally favored by determinations of the mixing angle dp from decay data. 
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